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SUMMARY 


Experimental and theoretical investigations of the effect 
of suspended solid particles on the performance of the compressor 
cascade is presented. The experimental investigation was carried 
out in a specially built cascade tunnel, using quartz sand 
particles. The cascades were made of NACA 65 ( 10)10 airfoils. 
Three cascades were tested, one accelerating cascade and two 
diffusing cascades. 

The theoretical analysis assumes inviscid and incompressible 
two dimensional flow. The momentum exchange between the fluid 
and the particle is accounted for by the interphase force terms 
in the fluid momentum equation. The modified fluid phase 
momentum equations and the continuity equation are reduced 
to the conventional stream function-vorticity formulation. 

The method treats the fluid phase in the Eulerian system and 
the particle phase in Lagrangian system. In addition, the 
inelastic collision between the particle and the blade surfaces 
are accounted for in the computation. The equations were 
solved numerically in a shifting coordinate system, which 
enables space marching solution of the vorticity equation. 

The experimental results indicate a small increase in 
the blade surface static pressures, while the theoretical 
results indicate a small decrease. The theoretical analysis, 
also predicts the loss in total pressure associated with the 
particulate flow through the cascade. 
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INTRODUCTION 


The ingestion of sand in aircraft engines can cause serious 
damage to the engine in terms of. sudden loss of power due to 
the momentum exchange between the fluid and the solid particles. 
This is only a temporary phenomenon and the effect lasts only 
for the duration of the sand ingestion. In addition to this 
temporary effect/ theJe is also a permanent irreversible 
change in the aerodynamic performance due to the erosion of 
various engine components. This effect will be discussed in 
another report. The temporary decrease in the power output 
is of serious concern to aircrafts operating in the desert 
environment where apart from the airborne sand ingested while 
on take-off and landing, the occasional sand storm and the 
resulting relatively high concentration of sand can pose a 
serious threat to the operation of the aircraft. Many such 
incidences have been reported. 

The phenomena of sand laden airflow through a turbomachine 
is treated as two phase flow. In most cases the solid particles 
mass concentration encountered is very low and subsequently 
the volume concentration is negligible. For example, even 
a mass of 2 gm/ff^ is equivalent to about six percent mass 

-3 

concentration and a volume concentration of 2.8 x 10 percent 
only. 

This low concentration enables one to neglect the particle- 
particle interaction and treat the individual particles. 

This approach was used by Tabakoff and Hussein [1] to compute 
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the particle trajectories and study the behavior of the 
individual parJ^.icles through a given turbomachine. In addition, 
they also included the inelastic collision between the particle 
and the solid surfaces in these computations by means of 
experimental correlations for the so-called restitution ratios. 
Experimental data for such- ratios on the impact of sand and 
ash particles for a variety of target metals are available [2, 3]. 
Tabakoff and Hussein’s approaches have led to successful 
computational models to predict the particle trajectories through 
turbomachines. They also studied the particulate flow influence 
on turbomachinery performance by using a simple one dimensional 
model [4] , It was assumed that two stream tubes exist in the 
flov; field around a given airfoil - one at the pressure side 
and the other at the suction side. The gas flow without 
particles was used to determine the nondimens ional area of 
these stream tubes and one dimensional gas particle (two phase) 
flow equations were solved along these stream tubes to obtain . 
the particulate flow pressure distribution along the airfoil. 

The above method was successfully applied by Tabakoff 
et al. [5] to an axial flow turbine. However, the particles 
used were corn cups which are lighter than sand particles and 
tend to follow closer to the fluid streamlines, while even a 
sand particle of 30-40 micron size does not follow the fluid 
streamline. The major limitation of this method stems from the 
fact that the assimption of particles following the fluid 
streamline is not valid for sand particles. 
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In order to fully undc<rstand the flow phenomena under 
sand laden conditions, experiments' were conducted on 
NACA-65 series airfoil profiles in a specially built 2-D 
cascade tunnel. In addition, a theoretical analysis of the 
two dimensional inviscid incompressible particulate flow 
through cascades was carried out, which gives a new insight 
into the various aspects associated with particulate flow 
field. 


DEVELOPMENT OF FLUID-PARTICLE FLOW EQUATIONS 

The governing equations for two dimensional gas particulate 

flow are derived under the following assumptions: 

(1) The particles are spherical in nature. 

(2) The volume concentration of the particles are negligible 
and hence there is no particle-particle interaction. 

(3) The particle-fluid coupling is only due to viscous drag 
forces and the coupling is one way, i.e. there is no 
momentum transfer from particle to fluid. The energy 
dissipation of the particle while deceleration is due to 
viscous wake and this momentum does not contribute to an 
increase in fluid total pressure. This is due to the 
fact that the mean free path of the particle is very 
large as compared -to the particle diameter (very low 
volume fraction, Soo, S.L., Ref. [6]). 
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Though the fluid particle interaction is due to viscous 
forces, the viscous effects are neglected in the fluid 
momentum equation and the fluid is tre-ited as inviscid. 

(5) The flow is incompressible and there is no heat transfer 
between the fluid and the particles. 

(6) The particles are treated in the Lagrangian system and 
the fluid in the Eulerian system. 


Eulerian Formulation of the Fluid Phase 

The continuity and momentum equations of the fluid phase 
in Cartesian coordinates may be written as follows : 
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where and are the interphase force term and represents 
the momentum transfer between the fluid and particles which 
will be treated later. 

Equations (2) and (3) can be reduced together with 
equation (1) to give: 
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where w is th^s vorticity defined as 
3u^ 3v^ 

5/-3^=“ . <"> 

The following relations for the definition of a stream function 
can be written; 


„ 3, _w ^ V = - li 

g t» Pg 3y ' g b Pg 3x 


( 6 ) 


where w is the mass flow through cascade, and b is the blade 
height of the cascade. Equation (6) satisfies the continuity 
equation (1) and equations (4) and (5) then reduce to the 
conventional stream function vorticity formulation. 

liil L. IJi + _1 x ( 

3y 3x 3x 3y ” p^k 3y p^k 3x ' 

3^^ + = iii ( 

3x^ 3y^ ^ 


where 



Such formulation of equations (7) and (8) are easier to handle 
in a body fitted coordinate system. One may define a coordinate 
transformation such that 


n = f^(x,y) , ^ = f 2 (x,y) 


(9) 


Then the above equations reduce to 
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where 


« = + y2 , 


J = x-Y - X y- , 
5^n 


3 = x-x + YrY / 
5 n ' 


Y = + Y^ 


(12) 


There are many choices for the functions f^ and ± 2 * 
Referring to Fig. 1, one can define that the lines of constant 5 
lines of constant x are the same, i.e. 5 = x. This will 
help to impose the periodicity boundary condition for the 
cascade flow, then 

x^ = 0.0 - , x^ = 1.0 , 

ct = / 3 = y^y^ , (13) 

Y = 1 + y^ / J = 


and the vorticity transport equation (10) is modified as 

Cii \p - 0 ) ip- = [-F + y F - y_F ] 

5 n n ? kp X n Yr ^ y 

g n -^5 n 

One can apply the transformation to the original momentum 
equations (2) and (3) and separate the two pressure terms, 
resulting in the following equations for pressure: 


( 14 ) 
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(16) 


From the solution of equations (11) and (14) # one obtains 
the flow field. Equations (15) and (16) are then used to 
evaluate the pressure by integrating these equations for the 
given inlet conditions to the cascade. The boundary conditions, 
the method of solution and the choice of the function 
n = f(x,y) are discussed in Appendix A. 


Lagrangian Formulation of the Particle Phase 

As mentioned earlier, it is assumed that only the fluid 
drag is the major force acting on the particles. The equation 
of motion of the single particle in Cartesian coordinates is 
given by: 


du 

m 

p dt 
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where and are the drag 

m is their mass, and u_, v_ 
P P P 


forces experienced by the particle , 
are particle velocities in x and 
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y directions. It is assiimed that the particles are spherical 
and the drag forces are given by: 
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+ Re) 


1.0 < Re < 4.0 


Cj^ = 21.9416 Re“°*^^® + 0.3240 4 < Re < 2000 

Cj^ = 0.4 2000 < Re £ 3 x lo'^ 


The Reynolds number in equations (19-22) is given by; 


( 22 ) 


^ ^ '"g ■ “p’ ^ <23 

Particle Dynamics 

The equations of motion for the particles are solved by 
time marching technique [7] . At each time step knowing the 
particle location, the flow properties are computed by inter- 
polation. Then the particles drag forces and accelerations are 
calculated. If the particle impacts on the solid boundary, 
then the exact impact point is located by Newton-Raphson 
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iteration technique. At this point the experimental correlation 
data for the restitution ratios is applied. 

The magnitude and direction of the particle rebounding 
velocity after collision with a solid surface must be known 
in order that the solution of the particle equations of motion 
be continued beyond the point of collision. An experimental 
study was conducted to investigate the particle impact- 
rebound characteristics. The following correlations were 
derived for 165 micron quartz particles impacting aluminum 
6064 metal surface. 

Normal restitution ratio 

= 0.993 - 1.76 6 + 1.56 6^ - 0.49 3? 

^n^ 1X1 

Tangential restitution ratio 
Vt 

= 0.998 - 1.66 + 2.11 3^ - 0.67 3^ 

^1 

where V and V are the impinging and rebounding particle 
velocity components normal to the impacting surface and 
V. and V. are the tangential components parallel to the 
impacting surface. The angle between the incoming particle and 
the solid surface is denoted by 3^. 

The solution of equations (24) and (25) at the point of 
impact gives the rebounding particle velocity components. The 


(24) 


(25) 
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solution of the particle equations of motion is continued 
from the impact point until the particles exit the downstream 
boundary of the cascade. 


Computation of the Interphase Force Terms 


The computation of the force terms at each grid point 
where the fluid solution is to be found is based on the method 
employed by Crowe [8] . The average interphase force terms 
can be evaluat'd at any grid point i,j by considering a 
cell ABCD surrounding the grid point i,j as shown in 
Fig. 2. 

Since the number of particles entering the cascade can 
• be very large, it is assumed that each trajectory line, I, 
computed through the cascade is associated with a particular 
mass fraction of X^. If the total particles mass entering 
the cascade is M^, then the total number of particles 
represented by any trajectory line is given by: 



M /m 
P P 


(26) 


where m^ is the mass of the single particle. 

The number density of the particles in cell i,j due to 
the trajectory I is given by: 



= n^ At/V ,, 
p cell. . 

f j 


( 27 ) 
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where volume, and Ax is the residence time of 

the particles inside the cell .ABCD- associated with the 
trajectory line I. 

The interphase force terms are 
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and are the average particle accelerations 

^if j 


associated with trajectory I through the cell i,j. 

are the sum of the 


The total force terms F and F 




forces associated with all trajectories at cell i,j, i.e.. 
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and 
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If : 


i 
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These forces are then used in the computations of the 
flow field with particle presence. 
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EXPERIMENTAL FACILITY 


The design of the experimental set-up to study particulate 
flow through cascades is complicated. The most important constraints 
are to insure a good distribution of the particle concentration 
and particle velocity on the entire test section. Such distribution 
of particle concentration is difficult to achieve in a horizontal 
tiinnel. In addition the constant area duct of the tunnel to the 
entrance of the test section needs to be of specified length to 
accelerate the particles to a reasonably high velocity. However 
such a long consteuit area duct is associated with a rapid 
boundary layer growth and hence is not well suited for the puirpose 
of cascade studies. Any attempt to provide boundary layer control 
with the use of porous tunnel walls is rendered futile , since 
the particles can block off the porous walls as well as erode 
the surface away. The problem was overcome by providing a 
sudden contraction just before the test section to remove any 
very low velocity fluid near the entrance to the test section. 

This sudden contraction proved to be of great help in providing 
a reasonably uniform flow at the cascade test section over 
about 80 percent of the test section area. 

A schematic configuration of the tunnel is shown in Fig. 3. 

It consists of the following components: particle feeder (A) , 

main air supply pipe (B) , settling chamber and particle injector 
(C) , accelerating tunnel (D ) , test section (E ) , and exhaust tank 
(F). The equipment functions as follows: a measured amount 
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of abrasive grit of a given constituency is placed into the 
paxrticle feeder (A; , The particles are fed into a secondary air 
source and blown up to the particle injector in the settling 
chamber (C) , where it mixes with the main air supply (B) . The 
particles are accelerated by the high velocity air in the constant 
area duct (D) , before impacting the cascades in the test section 
(E) . Past the test section the particulate flow is exhausted 
through exhaust collector (F) . 

Since the particles are accelerated in the constant area 
duct by the aerodynamic drag forces, their velocity before 
impacting the cascade would depend upon the air velocity, the 
particle size and the length of the acceleration section (D) . 
Figure 4 gives an illustration of the dynamics of relatively 
large 165 micron particles with air flow velocity of 130 m/s. 

From this figure it can be seen that the particles final 
velocity is an exponential function of the tunnel length. 

Based on these findings , a tunnel length of 3 meters was 
used in obtaining the experimental data. The test section 
(E) of rectangular cross section 37.5 mm x 135 mm, was slightly 
smaller in cross section than- the accelerating duct (D) and is 
fitted with transparent walls. High speed photographic technique 
was used to study the sand particles in this section. In 
this manner the velocity of the approaching sand particles was 
obtained and compared to the theoretical predictions. In 
addition the particle distribution was studied by sampling the 
flow in different locations. 
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Instrumentation 


The primary and secondary air mass flow rates were measured 
using standard ASME orifice flow meters. The blade surface 
pressure distribution, the total pressure in the settling chamber, 
and the test section wall static pressures were measured using a 
sccuiivalve-pressure transducer-digital recorder system. 

RESULTS AND DISCUSSIONS 
Experimental Results 

For the purpose of comparison experiments were conducted 
on three different cascades according to NACA 65- (10) 10 airfoil 
specifications [9]. The airfoils were made of 6064-T6 Aluminum 
and produced by extrusion process. The chord of these ar rf oris 
were 50.8 mm. These airfoils were tested at a space chord 
ratio of 0.5. The low space chord ratio was chosen to ensure a 
minimum of five blade passages in the tunnel cascade. All the 
measurements were taken in the cascade center passage. The 
locations of the surface pressure taps are as shown in Fig. 5 
and their axial positions are as given in Table I. 

The experiments were conducted on three different cascades. 
One was an accelerating cascade and the other two were diffusing 
cascades. The particles used were quartz sand of 165y mean 
diameter. The details of the cascade inlet and exit angles and 
the stages and incidences are as shown in Table II. However 
it was observed that the stagge"*: setting of the cascade II was 
wrong by 3® and the air inlet angle was wrong by 4®. The 
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air inlet angle changed, since the guides used to locate the 
cascade were eroded by the sand resulting in improper positioning 
of the cascades with reference to the air flow direction. This 
problem was later corrected in the testing of the third cascade. 

All the experiments were conducted at an inlet air velocity 
of 130 m/sec, which corresponded to a particle velocity of about 
100 m/sec. Figure 6 shows the surface pressure measurements in 
absolute quantities as obtained in a particular test. 

Inspection of this figure shows that the overall static pre- 
ssure level increases in the entire cascade with the injection 
of sand in the tunnel. This is from the increased resistance 
to the flow due to the presence of the particle in the piping 
system. In addition, attempts made to measure the total 

pressure with the particles flowing proved to be futile, since the 
particles clogged up the probes. The only pressure measurements 
that were carried out were the blade surface pressures in the 
cascade passage, the inlet and the exit wall static pressures. 
However it was noted that the inlet and exit v/all static pressures 
increased nearly by the same amount with the presence of the 
particles. Since the air mass flow through the tunnel did not 
change, and the change in the air density at the test section 
was negligible, the inlet total pressure was computed as the 
sum of the inlet dynamic head and the inlet wall static pressure. 

For a better understanding of the results , they are presented 
in the form of non dimensional pressure coefficient, S, which is 
defined as; 
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where P, is the inlet total pressure, 

^i 

Pg is the surface static pressure, 

P,. is the inlet static pressure, and 

Jm 

is the inlet dynamic head. 

Figures 7 through 9 show the observed results. From the inspection 
of these figures , it can be seen that the change due to particulate 
flow, in the surface pressure distribution measured does not 
follow a specific pattern. In the case of the accelerating 
cascade. Fig. 7, there is a uniform change on the blade suction 
side and there is a significant change in the rear half of the 
blade pressure side. These changes are in the same direction 
cind the decrease in the pressure coefficient indicates a higher 
surface static pressure. In the case of the diffusing cascade, 
it can be noted from Figs. 8 ?nd 9 that the overall tendancy 
is a decrease in the pressure coefficients on both surfaces . 

Since, for most part of the blade surfaces, the pressure coeffi- 
cient with particulate flow decreases , one can conclude that 
there is an increase in the surface static pressures. 

Theoretical Results 

Using the analytical approach described earlier, the results 
for the same three cascades are presented in Figures 10 through 23. 
The trajectories through accelerating cascade for three different 
size particles 5, 10 and 165 microns in diameter are shown in 
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Figures 10, 11 and 12 respectively. Figures 13, 14 and 15 show 
the trajectories of the 5, 10 and 165 micron particles through 
the diffusing cascade III, From these figures it Ccui be seen that 
the large particles are least influenced by the air, while the 
small particles are highly influenced. However, the particles 
are more influenced in the case of the accelerating cascade 
than in the case of the diffusing cascade. In addition, it can 
be observed that there is a region of no particles on the rear 
half of the blade suction surface. It can also be concluded 
from these figures that the particle local concentration can be 
much higher thsui the mean concentration, particularly near the 
blade pressure surface. 

Figures 16 and 17 show the streamline patterns with and 
without particles, for two cascades. In both cases, the con- 
centration of particles is 10% and the particle dieoneter is 
165 microns. The streamlines with particulate flow bend away 
from the blade pressure surface towards the suction surface. 
However, these changes can be clearly seen only close cind 
beyond the trailing edges of the cascades. The maximum change 
was only of the order of few percent. The movement of the stream- 
lines away from the blade pressure surface results in a decrease 
in velocities on the pressure surface and the movement of the 
streamlines closer to the suction surface results in an increase 
in the velocities on the suction surface. 

Figures 18, 19, 20 and 21 show the effect of the particles 
on the cascade pressure distribution. The change in the blade 
surface pressure coefficients for the three cascades considered. 
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at identical particle and air velocities for a =* 0.10 and 
dp =s 165 microns, are shewn in Figures 18, 19 and 20. The effect 
of the particles show decrease on the blade surface static 
pressures, especially on the blade suction surface. The pressure 
distribution on the blade pressure surface is more sensitive to 
the particle presence and does not follow an exact pattern in the 
case of these large particles. Figure 21 shows the effect of 

small partices, d = 10 microns, on one of the diffusing cascades. 

P 

It can be seen from this figure that for very small particles, 
the pressure coefficient increases on both surfaces indicating 
a decrease in the surface static pres’=5ure on both surfaces. 

A simple explanation for this could be that there are two different 
phenomena happening inside the blade passages. First, the 
increased concentration of the particles near the blade pressure 
surface, results in a considerable decrease in the total pressure 
and velocities of the fluid near the blade pressure surface, 
which ar^; accompanied by an increase in the fluid velocities 
on the blade suction surface. However, the fluid velocities 
are altered in such a way that the static pressure variation 
in the flow field satisfies the velocity distribution in the 
flow field. For this reason, the theoretical method predicts 
an overall decrease in the static pressure in the flow field. 

A comparison of the theoretical and experimental results 
indicates that there is small disagreement. The experimental 
results indicate small increase in the surface static pressures 
for the most part of the blade profile, while the theoretical 
results predict a small decrease in the surface static pressures. 
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There are many deficiencies both in the theoretical model 
cuid in the experimental techniques. The theoretical model, 
assumes incompressible flow, but the real flow is compressible. 

In addition, the theoretical model, assumes sharp leading and 
trailing edges. In the case of the experimental technique, the 
inlet total pressure could not be measured exactly and the inlet 
static pressure increases on introducing the particles. From 
the particle trajectories through the cascade one can observe 
an increased particles concentration near the blade pressure 
surfaces and it is possible that the local concentration can 
be suficiently high so that there is some contribution to the 
surface pressure changes. Furthermore, there is erosion taking 
place continuously as the testing progresses. In order to reduce 
any contribution from this aspect, the concentration has to be 
very low. However, at very low concentrations, the change in 
the pressure distributions can be so small, that it csui not be 
measured. In addition, the measured surface pressure distributions 
are time averaged values using the sc^nivalve system. There 
could be considerable interaction between the pressure measure- 
ment system and the continuously changing surface static 
pressures. The tests on casc&de III were carried out using 
multitube manometer instead of the scanivalve system. During 
these tests it was observed that the surface static pressure 
was continuously fluctuating during the entire testing period 
on introduction of the particles. All of the above mentioned 
problems could be the reason why the small difference exists 
between the experimental and theoretical results. 
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The mean exit total pressure of the cascade is always 
less than the mean inlet total pressure under particulate 
flow conditions . The theoretical analysis was carried out for 
various particle sizes from 10 to 250. microns diameter and 
particle concentrations of 5 to 20 percent. 

Based on the mean exit total pressure, a loss coefficient, 
C, due to the presence of particles can be defined as 
follows : 




particle 



where P = mean inlet total pressure, , 

= mean exit total pressure, 

= mean inlet dynamic head. 

Figure 22 shows the effect of particle diameter on the 
total pressure loss for of the cascades considered for the 

same particle concentration. The total pressure loss is much 
higher for the case of the accelerating cascade and the total 
pressure loss decreases rapidly as the particle size increases 
to a very low level. Under identical conditions, the fluid 
drag is proportional to the total surface area of the particles. 
For the same particle concentration, the total surface area is 
proportional to 1/^p* For this reason, the small particles 
give rise to a large decrease in total pressure. In the case 
of the decelerating cascade, the difference in velocities 
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between the rebounding particles and the fluid is not as high 
as in the case of the accelerating cascade. This leads to the 
loss coefficient in decelerating cascades being normally lower 
than that in accelerating cascades. The effect of concentration 
on the loss coefficient is shown in Fig. 23 for the two cascades 
considered. It can be observed that the loss coefficient increases 
linearly with the particle concentration. In order to enable 
comparison under identical conditions , the particle velocity to 
air velocity ratio at the inlet of the cascade was the scime in 
all the cases. In a real situation, the velocities of the 
smaller pairticles will be close to the air velocities. Consequently, 
the total pressure loss coefficients may not be as high as indicated 
in Pig. 22 for the small particles. 

CONCLUSIONS 

The experimental investigation on particulate flow through 
air foils in cascade shows a small decrease in blade surface 
pressure coefficients with the presence of particles. The 
theoretical analysis indicates that the streamlines bend towards 
the suction surface with the introduction of the particles. 

Though the theoretical analysis also shows only a small change 
in the pressure coefficients with the introduction of the particles , 
the pressure coefficients were found to increase with the presence 
of particles. The reason for -^his difference could be mainly 
from the interaction between the pressure measurement system 
and the effect of particles on the flow conditions in the cascade 
tunnel like increased turbulance level and back pressure. The 
theoretical analysis also predicts the total pressure loss 


associated with the particulate flow. This loss is found to be 
directly proportioned to the particle concentration. The total 
pressure loss is considerably higher for smaller particles than 
bigger particles \inder identical inlet conditions to the cascade. 
In a real situation the ratio of particle velocity to the air 
velocity will be higher for the smaller particles as compared 
to the big particles, and the total pressure loss may be slightly 
different than that used for comparison. However, this loss 
in total pressure can significantly alter the pressure ratios 
in a multistage machine, leading to a considerable decrease in 
the performance. 
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NOMENCLATURE 



D 

F 

(if j) 
I 
J 
k 

9 

in 

m 

P 

q 

Re 

S 

At 

u 


V 


V 

p 


w 

(x,y) 


blade height in cascade, (m) 
coefficient of drag 
chord, (m) 
axial chord, (m) 
diameter, (m) 
drag force, (N) 
interphase force, (N/m ) 
grid point indices 
index of the trajectory 
transformation parameter 
constant defined as W/bp 

g 

mass flow rate of particles along a trajectory (kg/sec) 
mass 

2 

pressure, (N/m ) 

12 2 
dynamic head, pV , (N/m ) 

Reynolds number 

pressure coefficient 

residence time, (sec) 

velocity in the x direction, (m/sec) 

average acceleration of the particle in the x direction, 
(m/sec^) 
volume, (m^) 

average acceleration of the particle in the y direction, 
(m/sec^) 

mass flow rate of gas through one cascade 
Cartesian coordinates passage, (kg/sec) 
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a transformation parameter or particle mass concentration 

S transformation parameter or angle 

Y transformation parameter’ 

C total pressure loss coefficient 

2 

li viscosity, (N sec/m ) 

(?,n) transformed coordinates 

ijj . stream function 
0 ) vorticity, (1/sec) 

Subscripts 
e exit 

g gas 

i inlet 

p particle 

t total 

(x,y) X and y components 

C / Ti partial derivatives with reference to 5 , D 
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TP iBLE r 

LOCATION OF SURFACE PRESSURE PROBES 


Suction 

Pressure 

x/C 
' ax 


*■1 

0.040 

^2 

^2 

0.125 

S 3 

^3 

0.250 

®4 

*“4 

0.350 

®5 

^5 

0.475 

®6 


0.625 

S 7 

^7 

0.725 

®8 

^3 

0.825 

S 9 

^9 

0.890 
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\ 

TABLE II 


DETAILS OF 

CASCADES USED 

IN TESTING 

PROGRAM 

Parameter 

Cascade I 

Cascade II 

Cascade III 

Air Inlet 
Angle, 3^ 

0.0 

•35° 

45° 

Stagger 

-20° 

+15° 

+25° 

Camber 9 

35° 

35° 

35° 

Incidence 

-3.00 

-3.00 

-3.00 

Aspect Ratio 

0.75 

0.75 

0.75 

Pitch-Chord 

Ratio 

0.5 

0.5 

'0.5 

No. of Blades 

6 

7 

9 
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FIG. 2. CONTROL VOLUME FOR COMPUTING FORCE TERMS 
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fig; 3, SCHERATIC OF CASCADE EROSION TUNNEL. 
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FIG. PARTICLE VELOCITY ALOIIG THE TUNNEL. 
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FIG. 6. MEASURED PRESSURE DISTRIBUTION OF CASCADE I. 
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FIG. 10. PARTICLE TRAJECTORIES THROUGH CASCADE I 
d„ = 5 MICRONS. 





E TRAJECTORIES THROUGH CASCADE I 
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FIG. 13. PARTICLE TRAJECTORIES THROUGH CASCADE III 
<5 - 5 MICRONS. 
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FIG. W. PARTICLE TRAJECTORIES THROUGH CASCADE III 
d„ = 10 MICRONS. 
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FIG. 18. COMPUTED PRESSURE DISTRIBUTION FOR CASCADE I 
FOR « = 0.0 AND a = O.'lO.' 
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FIG, 23. COMPUTED EFFECT OF PARTICLE CONCENTRATION 0! 
TOTAL PRESSURE LOSS COEFFICIENT. 
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APPENDIX A 


The method of solution of the stream function-vorticity 
equations depends mainly on the choice of the coordinate system. 

A proper choice of such system simplifies the application of 
boundary conditions. The ^ coordinates were chosen as lines 
of constant x, which enables one to apply the periodicity 
boundary conditions for cascade flow. The ri coordinates 
were chosen as any body fitted coordinate lines. 

In the present problem, the vorticity generation mechanism 
is due to the interphase force gradients. In order to get a 
fairly smooth force gradient at all points, the grid system has 
to be fairly coarse. Since the error in the computed vorticity 
at the body surface can be large in a coarse mesh, it is 
preferrable not to use th^ vorticity at the body surface during 
the computation. This can be achieved by exploiting the 
transport properties of the vorticity equation (10) . In the 
case that the lines of constant n are the same as the streamlines 
the vorticity equation can be solved by marching along the 
streamlines. Since the gradient ijj^ is zero, the term w ip- 
vanishes. However, one does not know the streamlines apriori. 

The other alternative is to neglect the term in 

equation (10) and march along lines of constant n* The error 
in the computed vorticity can be minimized if the n coordinates 
are updated to the new streamline pattern at every iteration. 

The advantages of this method are that the only information 
needed is the vorticity values at the inlet boundary and the 
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solution procedure does not use the vorticity at the blade 
surfaces. 


Neglecting the term on the left hand side of 

equation (10) results in 






(A.l) 


The stream function equation can be simplified to 

2 




(A. 2) 


where 


A = 


3y 


3C 


T1 


3ti 


B = 


3y„(i + Yr) 3y 


3n 


C 

35 ' 


(A. 3) 


C = y 


0 ' 


D = (1 - y^) f E = -2. * y^ 


Boundary Conditions 

The boundary conditions for the above stream function 
equation (A. 2) are (Fig. 1); 

For both the upstream boundary (A-B) and downstream boundary 
(E-F) , the inlet and exit angles and ^2 specified as: 

V 


tan 3i = (^) = (Yr-Y 

^ g inlet ^ 0^0 ^^let 


(A.4) 


tan = i:^) 

g exit 




(A. 5) 


exit 
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Bg[[ 


For the periodic boundaries (AH) , (BC) , (AF) and (DE) : 


^BC) = ^(AH) + 1-0 

(A. 6) 

^(DE) “ ^(GF) +1-0 

Along the airfoil surfaces: 

Suction Side (H-G) 

\jj = 0.0 (A. 7) 

Pressure Side (C-D) 

ip = 1.0 (A. 8) 

Trailing Edge Condition 

During the marching solution of the vorticity equation (A.l), 
one needs the information at the trailing edge vorticity. In 
potential flow, the vorticity at the trailing edge is zero. 
However, for the present case, where the total pressure is not 
uniform in the flow field, the trailing edge vorticity cannot 
be zero. The usual closure conditions applied for potential 
flow with uniform total pressure are that all the flow pro- 
perties are continuous and subsequently, the velocity gradient 
normal to the trailing edge is zero, or that the pressure side 
and suction side velocities are equal at the trailing edge. In 
the present case, where the total pressure is not the same on 
either blade side, the flow requirements are that the mass flow 
and static pressure are continuous across the blade trailing 
edge. This implies that the velocity field across the trailing 
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edge is not continuous, consequently, the vorticity is not 
zero at the trailing edge. The vorticity at the trailing 
edge was computed as the sum of the blade suction and pressure 
side vorticities [9] . 

Finite Difference Scheme 

Stream function equation ; 

The finite difference form of the stream function equation 

(A. 2) was obtained using the nine point central differences of 
2 2 

accuracy 0(An ,A^ ). Referring to Fig. (A.l) for any point i,j: 


A.i];. , .,+A-iI;. . .+ A- . .,_+A. il;. .t+Ac-tL. 

1 ^1-1, 3-1 2 ^1-1,3 3 ^1,3+1 4 ^i,3-l 5 ^ 1,3 


^6 ^i,j+l ^7 '^i+1,3-1 ^8 ^i+l,j ^9 '^i+l,j+l 


J . .03 . 

1/3 ifj 

^i/ j 


(A. 9) 


where 


= 


T.’ 


(4AnA^AQ) ' 


^2 (2A?) ' 


A^ = - A^ 


""4 “ UAn) ^ 


A^ = 1.0 


A = (— — — + — — T>r— ) /A 
^6 ^(2An) 


A^ = - A. 


^ ^ ^ ) a 

^8 ^( 2A ?) ^ ^0 


A 9 A^ 


and 


■^0 2 2^ 
^ AC An'^ 
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and the constants A, B, C, D and E are the transformation 
parameters defined by equation (A. 3) evaluated at point Ci/j). 

This gives raise to a system Of n equations where n is 
the number of .points at which \p is to be evaluated; 

[A] [^] = M (A. 10) 


The above system of equations can be solved by any known 
methods. For the present purpose, they were solved by a 
banded matrix solver (Refs, [10, 11]). 

Vorticity transport equation ; 

Considering equation '(A, 2), we can write it as 


CO . . — (0 . , 

i,j i~l,D ^ 
AC 


^ Prr 'P 


[-F, 


+ F 


g 

^ '2 t J 


“ F * Y tr 1 

^±-h,3 


(A. 11) 


which can be solved by marching along' the lines of constant n 
(i.e, the streamlines) 


to. . = to. 




[-F. 


i,j i-l,j k p u. , . ^ -xTi 

g 1 '2/j 


+ F r * y 

i-h,3 


- F * y,. 


(A, 12) 
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Solution Procedure 


The solution was started with an assumed pattern of 
streamlines satisfying the inlet and outlet flow angle require- 
ments. The vorticity was assumed to be zero and the irrotational 
flow streamlines were found for 21' streamlines. This was 
followed by the computation of trajectories for 100 particles 
introduced at the inlet. The particles equation of motion was 
solved by time marching technique. The particle impact point 
was found by Newton-Raphson iteration method. This was 
followed by the computation of the force terms, which were 
assumed to be associated with the streamlines rather than the 
fixed spatial coordinate. 

Once the force terms are computed, the vorticity is computed 
by marching along the n coordinates. This is followed by the 
stream function solution. Then the coordinates are updated 
and the process is repeated until the streamlines and the 
vorticity converge. This criteria was 0.01% relative error 
for the vorticity and 0.01% of the pitch for the streamlines. 
However, additional iterations were needed starting from the 
computation of the trajectories, since the force terms were 
associated with the streamlines rather than a fixed spatial 
point. The Flow Chart for the method of solution is given in 
Fig. A. 2. 

The pressure solution was obtained by integration of the 
momentum equation along x] and C directions . 
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FIG. A.l. FINITE-DIFFERENCE NOMENCLATURE. 
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FIG. A2. FLOW CHART OF THE COMPUTER PROGHAM. 
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